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PACS 05 . 40 . Fb - Random walks and Levy flights 

PACS 05. 10. Gg - Stochastic analysis methods (Fokker-Planck, Langevin, etc.) 
PACS 02 . 50 . Ey - Stochastic processes 

Abstract. - Phase transitions and effects of external noise on many body systems are one of 
the main topics in physics. In mean field coupled nonlinear dynamical stochastic systems driven 
by Brownian noise, various types of phase transitions including nonequilibrium ones may appear. 
A Brownian motion is a special case of Levy motion and the stochastic process based on the 
latter is an alternative choice for studying cooperative phenomena in various fields. Recently, 
fractional Fokker-Planck equations associated with Levy noise have attracted much attention 
and behaviors of systems with double-well potential subjected to Levy noise have been studied 
intensively. However, most of such studies have resorted to numerical computation. We construct 
an analytically solvable model to study the occurrence of phase transitions driven by Levy stable 
noise. 
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Introduction. — Phase transitions in coupled non- 
linear dynamical systems are one of the main topics in 
physics, especially in statistical physics, and have at- 
tracted much attention not only in natural sciences but 
in social sciences, since such studies are required to un- 
derstand various types of cooperative phenomena in en- 
sembles consisting of elements interacting with each other. 
Especially for systems with mean field couplings, owing to 
the simplification of interactions, phase transitions includ- 
ing nonequilibrium types of phase transitions that occur 
as a consequence of ergodicity breaking are investigated 
in detail and bifurcations of temporal attractors such as 
limit cycle and chaos with changes in the noise amplitudes 
have been found [1-3] (see also [4]). The starting point of 
such studies is usually a set of Langevin equations driven 
by white Gaussian noise, more precisely Brownian noise. 
By taking advantage of the self-averaging properties for 
mean field couplings, the master equation describing the 
temporal evolution of the empirical probability density of 
the system becomes a nonlinear Fokker-Planck equation 
(NFPE) [1,2,5-8] in the thermodynamic limit. The origin 
of the nonlinearity in NFPE with respect to the proba- 
bility density is the mean field coupling term which gives 
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rise to averaging over the empirical probability density in 
the thermodynamic limit. Especially for the nonlinearly 
coupled stochastic dynamics with quadratic potential, the 
form of the NFPE with the self-consistently determined 
mean field coupling term coincides formally with the form 
of the linear Fokker-Planck equation corresponding to the 
well-known Ornstein-Uhlenbeck (OU) process, i.e., the OU 
process of onebody system. Since the probability density 
for the OU process after large times takes a Gaussian form 
according to an H-theorem, the temporal evolution of the 
order parameters is analytically and easily obtained for 
the mean field coupled system with quadratic potential. 

Levy stable noise [9] is considered to be as important as 
Brownian noise in dynamical stochastic systems. Brown- 
ian noise is a special case of Levy stable noise, and Levy 
noise distribution exhibits the asymptotic power-law decay 
whereas Brownian noise distribution decays exponentially. 
Recently anomalous diffusion processes associated with 
Levy stable noise have attracted much attention in a va- 
riety of fields not only of natural sciences such as physics, 
biology, earth science, etc., but of social sciences such as 
risk management, finance, etc. In the context of physics, 
the microscopic description of the Levy processes associ- 
ated with anomalous diffusion has been proposed in [10]. 
Levy flights have been realized by the displacement of 
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the random walker subjected to Gaussian multiplicative 
noise as a random force. In such a system, the statistical 
property of the displacement of the random walker de- 
pends on its " memory" at small time scale, while it obeys 
fully Markoffian process, which reproduces a Levy flight, 
at large time scale. Furthermore, one of interesting appli- 
cations in physics can be found in subrecoil laser cooling of 
atoms [11], and the studies on this system show the valid- 
ity of use of Levy statistics in physics. To study subrecoil 
laser cooling process, the quantum mechanical analysis on 
the interactions between atoms and photons, which is sig- 
nificantly difficult in the system of dimension d > 2, is 
required. On the other hand, dealing with this problem as 
a classical inhomogeneous random walk (anomalous dif- 
fusion process owing to the long lifetimes of atoms with 
constant momentum close to zero) of 'dressed atoms' in 
momentum space allows us to use Levy statistics to easily 
obtain the behaviors of cooled atoms. This approach is 
shown to be equivalent to more complicated microscopic 
quantum approach, and a good agreement between theo- 
retical and experimental results is reported. Furthermore, 
Levy noise was studied from the viewpoint of improving 
neural signal detection [12]. 

Levy statistics, in general, is expected to be a powerful 
tool to investigate dynamical stochastic processes driven 
by long tail distributed noise, the superposition of a suffi- 
ciently large number of which does not follow the central 
limit theorem. Such a superposition is known to obey a 
Levy stable distribution. This fact is regarded as the gen- 
eralization of the central limit theorem. 

While the stationary probability density does not take a 
Gibbsian form in the system with Levy noise, it is shown 
that an H-theorem applicable to systems with Levy noise 
holds [13]. Surprisingly, the H-functional for such systems 
also takes the form of the relative entropy, the same form 
as the one for the systems with Brownian noise. Consider- 
ing such correspondences in spite of the significantly differ- 
ent statistical behaviors between Brownian noise and Levy 
noise, it will be important to investigate how statistical 
properties, especially cooperative phenomena in systems 
with Brownian noise is generalized and extended in sys- 
tems with Levy noise. The main target of this paper is to 
elucidate how behaviors of the nonlinear coupled stochas- 
tic dynamics with mean field couplings in the presence of 
Levy noise differ from those in the presence of Brownian 
noise, in terms of the existence of the equilibrium phase 
transition, stability condition for equilibrium solution, and 
the criticality. The phrases " bifurcation" or " phase transi- 
tion" have been used to describe the statistical properties 
of onebody stochastic systems in some literatures [14-16]. 
The "bifurcations" associated with Levy noise have also 
been investigated to understand the difference of the influ- 
ence on systems between Levy and Gaussian noise [15,16]. 
There, however, the "bifurcations" only imply the change 
of the form of the (stationary) probability density and 
are not bifurcations in the mathematical sense. The bifur- 
cations presented in this paper are totally different from 



them. We want to be concerned with genuine phase transi- 
tions or bifurcations where the stability exchanges between 
attractors occur. We will present rigorous analyses to ob- 
tain exact results in this paper. 

Model. — We deal with the following stochastic dy- 
namics of N one-dimensional elements: 

dxj(t) = m j ({x i (t)},t)dt + a({x i (t)},t)dL j , (1) 



where m and a 
and time t, Lj 



denote functions of x\{t), ■ ■ ■ , XN(t) 
J} is a Levy stable motion and dLj 
is determined by its characteristic function involving 
four parameters (a, (3, 7, D): E [exp (ikdLj)] = 
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o(di) with 



u)(k, a) — tan ^ for a ^ 1, u{k, a) = — - In \k\ for a = 1. 
E[-] denotes the average over the Levy noise dLj. Here 
the parameter a (0 < a < 2) characterizes the asymp- 
totic tail of the Levy stable distribution f(L; a, /?, 7, D) 
for a < 2 as f(L; a, (3, 7, D) ~ \L\-"- 1 with \L\ > 1. For 
a = 2, the Levy stable distribution reproduces a Gaus- 
sian distribution. The parameter (3 (— 1 < (3 < 1) is the 
skewness parameter defining the degree of asymmetry of 
the stable distribution, 7 (—00 < 7 < 00) is the center or 
location parameter which denotes the mean value of the 
distribution when a > 1 and D (0 < D < 00) is called 
the scale parameter which represents the generalized dif- 
fusion coefficient. The function to(k, a) is discontinuous at 
a = 1. The alternative choice for u(k, a) continuous with 
respect to a has been presented in some literatures [17,18]. 
However we will use the standard choice for uj(k, a) in this 
paper. 

In this paper we consider that the Levy stable noise 
dL, (i)'s are white noise and statistically independent with 
respect to the site j, i.e., they follow i.i.d. 

For the purpose of dealing with an analytically solvable 
model, hereafter the parameters a, [3, 7 and D are as- 
sumed to take same value for every site j and we treat the 
case of 
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mj (x(t),t) = -axj{t) + ^^F(s<(i)) 



(2) 



with a positive constant, J constant, F(x) an arbitrary 
bounded nonlinear function and a in eq. (1) a positive 
constant. Eq. (2) is well known to describe analog neu- 
ral network models with F representing a transfer func- 
tion, which is usually taken to be a bounded function 
such as a sigmoid function. A choice of non-monotonic 
functions for F was shown to improve the performance 
of associative memory neural networks [19]. It is also 
noted that the linearlity of mj in the variable Xj makes 
the model analytically tractable [2] as shown below. In 
this model, since the self-averaging property or the law of 
large numbers will hold for the mean field coupling term 

M ( t ) = W^2j F { x j{t)) in tnc lar g e N limit : the dynam- 
ics of the A^-body system (1) is expected to be reduced to 
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dealing with the onebody stochastic dynamics: 

dx(t) = m(x(t),t)dt + adL, 
m(x(t),t) = -ax(t) + JM(t) . 



(3) 



Notice that the dynamics (3) in the case of a = 2, i.e., 
in the special case where the stochastic process L is a 
Brownian motion, becomes the well-known OU process. 
This type of onebody stochastic dynamics driven by Levy 
stable noise (3) without the mean field interaction term 
M has been investigated in [20]. The important differ- 
ence between the model presented in [20] and our model 
is the presence of the coupling term M in our model, which 
makes the time evolution equation of the probability den- 
sity p(x, t) nonlinear owing to M(t) — F(x)p(x, t)dx. 
Hence phase transitions may appear in our model. 

Fractional Fokker-Planck equation and station- 
ary solutions. — The FPE corresponding to the 
stochastic differential equation with Levy noise takes the 
form of fractional FPE [13,21-23]. According to [22,23], 
the nonlinear fractional FPE (NFFPE) associated with 
the stochastic dynamics (3) will be given as 
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— [m(x, t) + 7ct] p(x, t\x ,t ) 
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where ^ r_1 (-) denotes the inverse Fourier transform and * 
represents the convolution. Using the Fourier transform of 
the probability density p(k,t) = f^° oo p(x, t) exp(ikx) dx, 
the FPE (4) is rewritten as 



% = -ak^+ik(JM + W )p 



- Da a \k\ c 



1 - ^]^| w ( fccr ' a ) 



p. 



(5) 



In the case of a =/= 1, one can straightfor- 
wardly show that the temporal solution for the 
FPE (5) after large times takes the form p(k,t) = 

and the pa- 



exp ^7efffc - D eB \k\ a (l - i/3 e ffjfj- tan 
rameters /3 e fj, 7 ff and D e g obey the dynamics: 



-aj eS + JM+ja, ^ = -aaD eB + Da a , £ {p cS D cH ) = 
—actD c f[l3 c f[ + D(3a a . In principle, M(t) is evaluated in 
terms of /3 e gf, 7 g and D cB by the definition of M(t) 
and the above-mentioned temporal solution p(k,t). For 
some special bounded functions F(x), the concrete form 
of M(t) can be represented explicitly in terms of an an- 
alytic function of these parameters. In particular, for 
F(x) = since, assuming j3 = 7 = and the initial con- 
dition p(x, t = 0) = S(x — xo), the temporal solution (for 
any t) is given by solving eq. (5) as 

{Da a \k\ a 1 
aa ^ ~ cx P( _aai )] + ikB (t) f j 



B(t) = x exp(-ai) + [ ds [JM(s)] exp [-a(t - s)] . 

Jo 

(6) 

Then the time evolution of the order pa- 
rameter M(t) is given by solving M(t) = 
exp [-Da a (1 - e- aat ) /{aa)\ sin B(t) together with 



dB(t) 
dt 



+ aB(t) — J exp 



Da a (1 



aa 



sin B(t) (7) 



with 5(0) = x Q [Shiino M. and Ichiki A., 8 th Asia-Pacific 
Conference of Fundamental Problems of Opto and Micro- 
electronics (Kokushikan University, Tokyo) 2008]. Here- 
after we will confine ourselves only to the case of a 1 
for simplicity. Then the stationary solution of (5) is ob- 
tained as 



Peq(k) 



exp 



i(JM + 70-) 



k 



D\k a a a f „k 7ra\ 

— u 1 - i(3-r- tan — 

aa V y fc 2 J 



(8) 



In general, the moments related to the stationary solu- 
tion (8) diverge. However, the mean field interaction 
term M is finite since F(x) is bounded. It is worth 
to mention that the temporal solution p(k, t) after large 
times corresponds to Levy stable distribution. This is 
reminiscent of a similar situation described by the H- 
thcorem for the stochastic processes with Brownian noise 
presented in [1,2], which ensures the convergence to a 
temporal probability density of Gaussian form after large 
times. The proof of the convergence to the Levy sta- 
ble distribution in our model will be presented in de- 
tail elsewhere. Notice that the solution (8) reproduces 
the equilibrium probability density for the usual OU pro- 
cess when a = 2. Using (8), the stationary value Mq 
of the temporal order parameter M(t) is determined self- 
consistently as the solution of the order parameter equa- 

,. /—£> a a \ . \ JMn+~i<j 1 Da a B tan(-7ra/2) 1 • 

tion Mn = exp ( ) sin - — <n - L — — — m 

the case of F(x) = sinx. The solution of this order pa- 
rameter equation is plotted as a function of the gener- 
alized diffusion coefficient D in fig. 1. We can see that 
phase transitions occur and stability of the order parame- 
ter changes with a change in the parameter D. The value 
of D when the nontrivial stable solution Mo 7^ vanishes 
with a gradual increase in D is analytically evaluated as 
D c = (aa)ln(J/a) in the case of F(x) = sinx, a / 1, 
/3 = 7 = 0, (7=1. In general, the critical value D Cl which 
is the value of D when the stability of Mq alters, is evalu- 
ated by solving the equation J J p cq {x) [dF(x)/dx] dx = a 
with respect to D according to the stability analysis for 
the stationary solutions mentioned bellow. Furthermore, 
it is found that the stable solution M 7^ scales as 

1 /2 

Mo cx (D c — D) 1 near the critical point D < D c in ac- 
cordance with the mean-field universality class, by using 
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Fig. 1: Solution of the order parameter equation with F(x) = 
sinx, J = 15, a = 1, a = 1, a = 1.2, 7 = 0, /3 = 0. The 
stable solution is plotted in the solid line and the unstable is in 
the broken line. The critical value of the generalized diffusion 
coefficient is evaluated analytically as D c = (aa/a a ) ln(J/o) ~ 
3.25. 



the stationary solution (8) and expanding the definition of 
the order parameter Mq with respect to the difference in 
the generalized diffusion coefficient D c — D. Notice that 
this scaling law also holds in the Brownian noise case since 
this law is independent of the parameter a. 

Stability analysis. — In what follows, we investigate 
the stability of the stationary solution Mo. Considering a 
small deviation of the probability density from the station- 
ary density Sp(x,t) = p(x,t) — p eq (x), the Fourier trans- 
form of this deviation turns out to follow the FPE 



dSp 



= -ak^f- + ik(JM + ja)Sp + ikJSMp cq 
ok ' 

Da a \k\ a tan™) Jp + 0(«P),(9) 



where SM is a small deviation of the order parameter 
M related to the deviation of the probability density dp: 
SM(t) = J F(x)Sp(x, t)dx. Here we consider the eigen- 
function of the Fokker-Planck operator in eq. (9), i.e., 
Such a solution of the FPE (9) is 



Sp\(k,t) 



iJk 



1 



J -{F') 



6M(t = 0)exp (Xt)p cq {k) , 



-A/a 



(10) 



where SM(t = 0) the deviation of the order pa- 
rameter M from the stationary solution Mo at 
time t = and (F') = J ^±p eq (x)dx, R = 
JdxJ^F(x)\k\-^ a p cq (k)exp(-ikx). Since 5p\(k = 
0, t) = for normalization condition, 1 — J(F')/(a+X) = 
is required for the case A/a > 0. Since we assume a > 
for physical interest, the positive eigenvalue of the Fokker- 
Planck operator A is given as A = — a + J(F'). Hence the 



condition for stable stationary solution M is 

a-J(F')>0. (11) 

Notice here that the a-dependence of the stability con- 
dition (11) exists only in the quantity (F 1 ), and this 
result reproduces the stability condition in the Brown- 
ian noise case when a — 2. For the choice of F(x) = 
sinx, one finds R = M , (F') = \ [p cq (l) + p oq (-l)] 
and then obtains the stability condition: < a — 

Jexp (-^) COS \lMo±2Z + /^tan(W2) l Xhc gtable 

and unstable stationary solutions in this case are plotted 
in fig. 1. 

Summary. — In summary, we have investigated how 
the properties of the nonlinear coupled stochastic dynam- 
ics with mean field couplings in the presence of Brownian 
noise is modified in the presence of Levy noise. We have 
shown the occurrence of bifurcations for mean field type 
NFFPE associated with the stochastic process which de- 
scribes nonlinearly coupled elements subjected to indepen- 
dent identical Levy noise. All of the analyses have been 
conducted rigorously in the thermodynamic limit for sys- 
tems under the influence of Levy noise with index a^l. 
Similar analysis can be carried out also in the case of 
a = 1. The detailed analysis including the case of a = 1 
will be reported elsewhere. In this paper, the temporal 
evolution of the order parameter has been obtained to- 
gether with the stability condition for its stationary value. 
We have also presented the scaling of the order parameter 
near the critical point. The exact results presented in this 
paper have been obtained by taking advantage of (i) the 
self-averaging property, or the law of large numbers owing 
to the presence of mean field couplings in the thermody- 
namic limit, (ii) the identical dynamics for every element 
constituting the system and (iii) the quadratic potential in 
the intrinsic dynamics. The choice of such potential gives 
rise to a linear system for which Levy distribution holds. 
The situation is reminiscent of the so called OU process 
with a Gaussian probability density that manifests itself 
in a linear system subjected to Brownian white noise. The 
case with a double- well potential will be of interest and is 
now under way. There also remains open problems of how 
fluctuation-dissipation theorems and critical fluctuations 
are modified under the condition that moments of Levy 
distribution in general do not exist. 

In literatures [15,16], the terminology "bifurcation" is 
used to describe the situation where the form of the proba- 
bility density of a single element subjected to a Levy noise 
changes between unimodal and bimodal functions. On the 
other hand, in this paper, genuine phase transitions or bi- 
furcations driven by Levy noise have been shown. Since 
in recent years the effects of Levy noise have been ex- 
tensively studied in various fields of science, a study of 
fundamental problems of phase transitions involving Levy 
noise should be of vital importance. Especially a rigorous 
analysis of such cooperative phenomena based on analyt- 
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ically tractable models is quite useful considering poten- 
tially wide applications of it. Detailed analyses including 
the derivation of eq. (7) and related issues as well as the 
case a = 1 will be reported elsewhere. 

* * * 
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